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Abstract 

It is argued why quasiparticle models can be useful to describe the ther¬ 
modynamics of hot QCD excluding, however, the case of a large number 
of flavors, for which exact results have been calculated by Moore. 

In a recent paper 0, Moore considered the thermodynamics of hot QCD with 
a large number of quark flavors, Nf 3> N c . Although this case is physically not 
directly relevant, it is still interesting to study the limit Nf —> oo as the next-to 
leading order result in the 1 /Nf expansion of the pressure can be calculated 
exactly as a function of the effective coupling G 2 = Nfg 2 / 2 in this model. 1 
Therefore, it is suggestive to use the large- Nf limit of QCD as a testing ground 
for the convergence of perturbative results 0 or resummation improved meth¬ 
ods 012 El From the strong-coupling behavior of the pressure at large Nf 
it was also concluded in 0 that the quasiparticle model 0 cannot be a good 
description for ‘real’ QCD near the QCD transition. In the following it is argued 
why such a conclusion has to be taken with some care. 

To start with, a sketchy derivation of the large- Nf pressure is given within 
the <I>-derivable approach. Setting the volume of the system to one, the QCD 
pressure for arbitrary Nf can be expressed in terms of the full quark and gluon 
propagators 0 0 (the ghost contribution and the subtraction of the vacuum 
pressure are not shown explicitly), 

p = Tr [In (—5 _1 ) + S S] - ^Tr [in (-D" 1 ) + UD\ + $[S, D\ . (1) 

The traces are taken over the 4-momentum and, accordingly, over the flavor, 
spin, color and Lorentz indices. 4> is the sum of all 2-particle irreducible bubble 
diagrams. The self-energies are obtained by cutting a corresponding line in 4>. 

1 Note that the resulting theory is not asymptotically free and leads to a Landau pole, so 

it is not defined in a strict sense. However, as argued in 0, the model can be understood as 
an effective theory with a controllable ambiguity which is small as long as the Landau pole is 
much larger than any other scale in the problem. 
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In the large -Nf limit, <L reduces to the 1-gluon exchange diagram for each flavor 
and contributes 4 Trim to p, which compensates the second term in the gluon 
contribution. While the gluon self-energy is of the order G 2 , the quark self¬ 
energy is of the order G 2 /Nf. Thus, the leading term in the 1 /Nf expansion of 
p is given by the free quark pressure, pi 0 = Trln(— S 0 1 ) ~ Nf. The term of the 
order Nf comes from the remaining gluon contribution [Tj, 

Pnio = -^Trln (-L> _1 ) • (2) 

Moore’s numerical computation for massless quarks shows that p n \ a first de¬ 
creases with increasing G 2 as expected, but then starts to increase again, ex¬ 
ceeding even the free value p n i 0 (G 2 = 0) before the coupling becomes eventually 
so large that the position of the Landau pole reaches the order of the temperature 
and the model becomes meaningless. 2 Clearly, such a non-monotonic behavior 
cannot be described by quasiparticles whose masses rise with the coupling. 

To understand the surprising behavior of the large -Nf pressure, and whether 
it might be relevant for real QCD, it is instructive to go back to an arbitrary 
number of flavors. Then, in dimensionless units and up to terms of order <? 4 , 
the perturbative expansion of the pressure reads |2j 

p( 3 \N c , Nf )= 7 -iL + l-5(N c + *N f ) f + (n c + \ Nf^j ' f , (3) 


with df = N c Nf , d a = N 2 — 1, and g = g/{ 47t). The normalization is such 
that the free gluon contribution is one. Note that to this order the result does 
not depend explicitly on the renormalization scale. Aside from the fact that 
in general the terms beyond 0(g 5 ) in the expansion of p cannot be calculated 
by perturbative methods, the series is not convergent. The best one can hope 
for is some sort of an asymptotic expansion. Such a series, however, should 
be truncated at a certain order, depending on the expansion parameter, to 
give the best approximation possible. Typically, one should truncate the series 
at the nth term when its modulus becomes as large as that of the preceding 
contribution, see, e. g., [EJ. Applying this prescription to the expansion © leads 
to the estimate that p^ is a reasonable approximation for g smaller than 


\/3 Ac + | Nf 
~L6 (N c + i N f ) 3G ' 


(4) 


p^ has a minimum at g* = | g% (which therefore might be physically relevant) 
with a relative depth of 


f _! P l3) g> 5 + (5) 

_ P^O) 144(4 + 7^/^) \k c +1n,J 

2 After submitting this article, Moore’s numerical results for p n \ Q as published in [3 were 
corrected in (7J- While the pressure still has a minimum, it now exceeds the free pressure only 
at a value of the coupling where it is already sensitive to the Landau pole ambiguity. 
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For g > t/ 3 , the 0(g 2 ) approximation, p^ 2 \ is expected to be more adequate 
until it obviously also becomes meaningless at 


(l 4 + 7 d f /d a \ 1/2 
V5 47V c + 5 N f ) 


( 6 ) 


The values of < 72 , 53 and 5 for N c = 3 and several numbers of Nf are given in 
table ??, and the representative case Nf = 0 is illustrated in figure^ 



O 

II 

i? 

£ 

II 

CO 

II 

Ci 

large Nf 

92 

3.2 

3.7 

3.9 

G 2 = 3.5 

93 

0.79 

0.96 

0.97 

G 3 = 3.4 

S 

0.0087 

0.0099 

0.0094 

0.14 


Table 1: The values of the coupling indicating the validity range of the quadratic 
approximation p^ of the pressure for several Nf and of p n i Q in the large- Nf 
limit, respectively, and the depth S of the minimum for the cubic approximation. 


Due to the estimated large validity range of the approximation p( 2 ' 1 one can 
expect that the monotonically decreasing behavior at large coupling is qualita¬ 
tively correct for the exact pressure. This is indeed what is observed in QCD 
lattice calculations with a small number of flavors 0, where the pressure de¬ 
creases monotonically with the temperature and reaches a rather small value 
at the transition temperature T c . Although it is not relevant for the large cou¬ 
pling which is of primary interest here, it is mentioned that the minimum of the 
approximation p^ is very shallow, and that the 0(g 4 ) contribution to p is for 
reasonable choices of the renormalization scale always smaller than the 0(g 3 ) 
contribution. 3 The latter observation might indicate that the perturbative series 
is not an asymptotic expansion in the strict sense. 

Considering now the large- Nf limit, one can read off the perturbative ex¬ 
pansion of pnio from eq. m after subtracting the term ~ TV/; with G = G/{ 47t), 




( 7 ) 


The perturbative structure of p n \ Q is rather different from the case discussed 
above. Although the (D{G 2 ) approximation vanishes at a similar value of the 
(effective) coupling, the estimated (as before) validity range of p[^ 0 is tiny, see 
table ?? and figure |21 Therefore, one cannot expect p n i 0 to become small at 
large coupling — contrary to the pressure at small Nf. Again, this expectation 
is in line with the exact result |T]. Numerically, the expression p^ o turns out to 

be a reasonable approximation 4 for values of G as large as the coupling at the 

( 2 ) 

minimum of p n i Q , while p K nl J Q never seems to be the optimal approximation. 

3 In figure^ the lower bound of the validity interval of p ( 3 ) is estimated from the magnitude 
of the 4th and 5th order contribution to p. 

4 It is noted that the G(G 3 ) contribution is (for relevant values of the coupling) always 

larger than the 0(G 4 ) contribution which is larger than the contribution of the order G 5 . 
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Figure 1: The perturbative approxi¬ 
mations, p and p( 3 \ of the pressure 
for Nf = 0, in units of the free value. 



G 

( 2 ) 

Figure 2: The approximations p nlo 
(3) 

and p nlo of the large- Nf pressure p n io, 
scaled by the free value. 


The truncation prescription for asymptotic series implies that the terms 
of a higher than a certain order, which depends on the size of the expansion 
parameter, almost cancel each other. For the pressure of QCD with a few flavors, 
in the region of physical interest, near the QCD transition, this occurs already 
after the leading order correction in g 2 . Hence, a quasiparticle model like jH], 
which incorporates the leading perturbative correction and resums higher order 
terms only partly to ensure thermodynamic consistency, appears to be justified. 
The situation is different for the large- Nf limit of QCD where the pressure p n i 0 
has a peculiar perturbative structure. This fact is physically plausible: Here the 
cubic plasmon term, i. e., the correction to the leading exchange contribution 
due to screening, is small since quarks screen less than gluons. Consequently, 
this term is relevant at much larger coupling than in the physical case. 

In conclusion, Moore’s argument — the quasiparticle model (j cannot de¬ 
scribe the minimum of p n i 0 at large Nf and is thus not expected to be a good 
description of the physical case either - would only be convincing if such a 
pronounced minimum also occurred in real QCD; otherwise the large -Nf limit 
is governed by different physics. By the argument that in both cases the strong¬ 
coupling behavior of the pressure is reflected in the appropriately interpreted 
perturbative results, or directly by contrasting the numerical results (lattice 
data for real QCD), this, however, seems indeed to be the case. 
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